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We present a spatio-temporal analysis of a statistically stationary rotating turbulence experiment, 
aiming to extract a signature of inertial waves, and to determine the scales and frequencies at which 
they can be detected. The analysis uses two-point spatial correlations of the temporal Fourier 
transform of velocity fields obtained from time-resolved stereoscopic particle image velocimetry 
measurements in the rotating frame. We quantify the degree of anisotropy of turbulence as a function 
of frequency and spatial scale. We show that this space-time-dependent anisotropy is well described 
by the dispersion relation of linear inertial waves at large scale, while smaller scales are dominated 
by the sweeping of the waves by fluid motion at larger scales. This sweeping effect is mostly due 
to the low-frequency quasi-two-dimensional component of the turbulent flow, a prominent feature 
of our experiment which is not accounted for by wave turbulence theory. These results question 
the relevance of this theory for rotating turbulence at the moderate Rossby numbers accessible in 
laboratory experiments, which are relevant to most geophysical and astrophysical flows. 


I. INTRODUCTION 


The energy content of turbulence is usually charac¬ 
terized by the energy distribution among spatial scales, 
either in physical or in Fourier space. For rotating, strat¬ 
ified, or magnetohydrodynamic turbulence [ij , waves can 
propagate and coexist with “classical” eddies and coher¬ 
ent structures, which advocates for a spatio-temporal de¬ 
scription of such flows. While temporal fluctuations are 
usually slaved to the spatial ones via sweeping effects in 
classical turbulence they are expected to be gov¬ 

erned by the dispersion relation of the waves for time 
scales much smaller than the eddy turnover time. The 
latter regime is the subject of wave turbulence theory, 
in which the assumption of weak nonlinear coupling be¬ 
tween waves allows to predict scaling laws for the spatial 
energy spectrum 0,0. 

It is a matter of debate whether wave turbulence the¬ 
ory (also known as weak turbulence theory) is a good 
candidate to describe rotating turbulence in the rapidly 
rotating limit. Solutions to the linearized rotating Euler 
equation can be decomposed into inertial waves, which 
satisfy the anisotropic dispersion relation 


.(k) = 2 oM, 


( 1 ) 


where Yl is the rotation rate and fc|| the component of 
the wave vector k along the rotation axis (referred to as 
the vertical axis by convention) 0. Accordingly, only 
fluid motions at frequencies a smaller than the Corio¬ 
lis frequency correspond to wave propagation. Fluid 
motions of weak amplitude and slowly varying in time 
{a 2n) can be described in terms of waves with nearly 
horizontal wave vectors: they tend to be two-dimensional 
and three-component (2D3C), invariant along the rota¬ 
tion axis, a result known as the Taylor-Proudman theo¬ 
rem. 


The trend towards two-dimensionality is a landmark in 
rotating turbulence, observed both in experiments and 
numerical simulations 01MI1- It originates from the 
modification of the nonlinear interactions by the Coriolis 
force, which yields preferential energy transfers towards 
modes with almost horizontal wave vectors [l2l - [l0 . In 
the frequency domain, this corresponds to the genera¬ 
tion of slow dynamics compared to the characteristic fre¬ 
quency at which energy is supplied to the system. These 
anisotropic energy transfers can be accounted for in terms 
of resonant and near-resonant triadic interactions of in¬ 
ertial waves [lfil - n0 . A major feature of rotating turbu¬ 
lence is the emer gen ce of inverse energy transfers in the 
horizontal plane |I3l [l9l - [^ . Inverse transfers between 
3D fast “wave” modes and the 2D slow “vortex” mode, 
mediated by near-resonant triadic interactions, are al¬ 
lowed at finite Rossby number only [H, [2^ . The 

2D mode is therefore fed either from the coupling with 
the 3D modes at finite Rossby number, or from direct 
energy input by the forcing. One naturally expects the 
energy within this 2D mode to undergo an inverse en¬ 
ergy cascade, similar to that of classical (non-rotating) 
2D turbulence [13, HH . 

Such coexistence between 2D and 3D flows is relevant 
to most experiments and numerical simulations, and can¬ 
not be accounted for by wave turbulence theory, which 
describes the direct energy cascade arising from resonant 
triadic interactions of 3D wave modes only [H, H^l ■ This 
theory therefore provides only a partial description of ro¬ 
tating turbulence in realistic systems, and careful exper¬ 
imental and numerical studies remain necessary to assess 
its range of validity. 

Laboratory experiments differ from most numerical 
and theoretical studies by the presence of rigid horizon¬ 
tal boundaries, where the rotatin g flo w achieves no-slip 
conditions through Ekman layers [SJ]. In a laminar Ek- 
man layer, the balance between the viscous and Coriolis 
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FIG. 1. (Color online) Experimental setup. An arena of 10 
pairs of flaps forces a turbulent flow in the central region 
of a water tank mounted on a rotating turntable. A laser 
sheet illuminates a vertical slice through a horizontal glass lid 
covering the fluid. 2D3C (two-dimensional three-component) 
velocity measurements are performed using stereoscopic PIV 
(particle image velocimetry) in a vertical square domain of 
size Ax X Az = 14 x 14 cm^, shown as a dashed square in 
panel (a). 


forces leads to a boundary layer thickness Seu — 

The belief is that, provided the experimental tank is deep 
enough, the bulk turbulent flow away from the top and 
bottom boundaries should resemble the one obtained in 
the ideal 3D periodic or stress-free domains considered 
in most numerical and theoretical studies. Closer to the 
horizontal walls, the boundary layers induce Ekman fric¬ 
tion that is not taken into account by most numerical 
studies. 

In the laboratory, the energy dissipation of rotating 
turbulent flows originates from three main contributions: 
bulk viscous dissipation of 3D flow structures, bulk vis¬ 
cous dissipation of quasi-2D flow structures (somewhat 
similar to the bulk energy dissipation of 2D turbulence), 
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FIG. 2. (Golor online) (a) Energy and (b) energy dissipation 
rate per unit mass for the 2D and 3D modes as a function 
of the rotation rate D. In both figures, the first data points 
(shown with arrows, at arbitrary abscissae) correspond to the 
non-rotating case D = 0. 


and dissipation through Ekman friction on the horizontal 
boundaries. 

In spite of the importance of the 2D mode in most 
geophysical and laboratory flows, the 3D fluctuations still 
play a crucial role in the dynamics of rotating turbulence, 
because they are much more efficient at dissipating en¬ 
ergy. This key feature is illustrated in Fig. [2 from data 
obtained in the present experiment (setup sketched in 
Fig. HI see Sec. m for details): we decompose the tur¬ 
bulent velocity field into a vertically-averaged 2D flow 
and a vertically-dependent 3D remainder, and show the 
corresponding energies and energy dissipation rates as a 
function of global rotation. For maximum rotation, al¬ 
though the 3D component contains a small fraction of 
the total kinetic energy, its dissipation rate is as large as 
that of the vertically-invariant 2D component. Moreover, 
both of these dissipations are larger than an estimate of 
the frictional losses due to laminar Ekman layers (see 
Sec. nni). An accurate description of the 3D structures 
of the flow is therefore essential to characterize the en¬ 
ergy fluxes in rotating turbulence at moderate Rossby 
number. 

A primary goal in this direction is to determine the 
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range of scales and frequencies for which 3D fluctuations 
follow the inertial-wave dispersion relation. This requires 
a full spatio-temporal analysis, which is very demanding 
in general for wave-turbulence systems: the accessible 
range of scales is usually limited in experiments, whereas 
long integration times are prohibitive in numerical sim¬ 
ulations. The case of rotating turbulence is particularly 
delicate because of the specific form of the dispersion 
relation O: the frequency is not related to the wave 
number, as in conventional isotropic wav e sy stems such 
as surface waves [s^ or elastic waves [H, [13, but to the 
wave vector orientation only. 


The recent studies of Clark di Leoni et al. [13 and 
Yarom and Sharon [13 constitute important steps for¬ 
ward in this respect. Using numerical simulation of ro- 
tating turbulence forced at large scale, Clark di Leoni et 
al. [38l | observe a clear concentration of energy along the 
dispersion relation of inertial waves, and provide a de¬ 
tailed analysis of the various time scales of the system. 
They observe a wave-dominated regime at large scale 
and a sweeping-dominated regime at small scale (see also 
Ref. El). In the experiment of Yarom and Sharon (39j | 
the forcing consists of a random set of sources and sinks 
at the bottom of a rotating water tank. They measure 
three-dimensional two-component (3D2C) velocity fields 
using a scanning particle image velocimetry (PIV) tech¬ 
nique, and observe also a good agreement with the iner¬ 
tial wave dispersion relation. In both Refs. [13 and [13, 
the inertial waves are observed at scales smaller than the 
injection scale, suggesting that they are fed by forward 
energy transfers, which is consistent with the predictions 
of wave turbulence theory. 


The aim of the present paper is to further analyze ex¬ 
perimentally the range of spatio-temporal scales at which 
inertial waves can be detected in rotating turbulence. 
Stationary rotating turbulence is produced by a set of 
vortex dipole generators which continuously inject tur¬ 
bulent fluctuations towards the center of a rotating water 
tank where measurements are performed. We showed in 
Ref. [13 that this configuration generates a double en¬ 
ergy cascade at large rotation rate: an inverse cascade 
of horizontal energy and a direct cascade of vertical en¬ 
ergy, which behaves approximately as a passive scalar ad- 
vected by the horizontal flow. Here we perform a detailed 
spatio-temporal analysis using two-point spatial corre¬ 
lations of the temporal Fourier modes computed from 
time-resolved two-dimensional three-component (2D3C) 
velocity fields measured by stereoscopic PIV in a vertical 
plane. We observe that, at large scales and frequencies, 
the spatio-temporal anisotropy of the energy distribu¬ 
tion is well described by the dispersion relation of iner¬ 
tial waves, whereas smaller scales are dominated by the 
sweeping of the waves by the energetic large-scale flow. 


II. EXPERIMENTAL SETUP 

The experimental setup, sketched in Fig. [1] is similar 
to the one described in Refs. [13, [4l| and is only briefly 
described here. It consists of a 125 x 125 x 65 cm^ glass 
tank, filled with H = 50 cm of water and mounted on a 
2 m-diameter rotating platform which rotates at a rate 
D in the range 0.21 to 1.68 rad s“^ (2 to 16 rpm). Tur¬ 
bulence is produced in the rotating frame by a set of ten 
vertical vortex dipole generators organized as a circular 
arena of 85 cm diameter around the center of the water 
tank. This forcing device was initially designed to gen¬ 
erate turbulence in stratified fluids, and is described in 
detail in Refs. [H, [13. Each generator consists of a pair 
of vertical flaps, 60 cm high and Lf = 10 cm long, al¬ 
ternatively closing rapidly and opening slowly in a cyclic 
motion of period Tq = 2TTjao = 8.5 s. The closing stage 
is achieved with the flaps rotating at an angular velocity 
O'/ = 0.092 rad s“^, and a random phase shift is applied 
between the generators. 

In the laminar regime, a single pair of flaps generates 
vortex dipoles with core vorticity ujf. Additional PIV 
measurements close to a vortex dipole generator indicate 
that this core vorticity is governed by the vorticity in the 
viscous boundary layer of the flap, w/ ~ (jfLf/5, where 
cTfLf is the flap velocity and 5 the viscous boundary layer 
thickness. In the present experiment, the vortex dipoles 
are unstable, and the closing of the flaps therefore pro¬ 
duces small-scale 3D turbulent fluctuations that are ad- 
vected towards the center of the arena by the remaining 
large-scale dipolar structure. 

The turbulent Reynolds number, computed from the 
rms velocity and the horizontal integral scale, is about 
400 in the center of the flow, and the turbulent Rossby 
number covers the range 0.30 — 0.07 for D = 2 — 
16 rpm [ 23 . 

We measure the three components of the velocity field 
u = Ux^x + UyGy + UzGz (with Bz Oriented vertically, 
along the rotation axis) in a vertical square domain of 
size Ax X Az = 14 x 14 cm^ located at the center of 
the circular arena at mid-depth, using a stereoscopic PIV 
system [3, [13 embarked on the rotating platform. These 
2D3C (two-dimensional three-component) velocity fields 
are sampled on a grid of 80 x 80 vectors with a spatial 
resolution of 1.75 mm. Two acquisition sets are recorded 
for each rotation rate D: one set of 10 000 fields at 0.35 Hz 
and one set of 1000 fields at 1.5 Hz. The combination of 
these two time series results in a temporal spectral range 
of three decades. 


III. 2D VS. 3D FLOW COMPONENTS 

In the present experiment, energy is primarily injected 
in the 2D mode (vertically invariant), but the instabilities 
in the vicinity of the flaps rapidly feed 3D fluctuations 
which are advected in the central region. Energy trans¬ 
fers between the 2D and 3D flow components, which van- 
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ish in the weak turbulence limit {Ro —> 0), are allowed in 
our system because of the moderate value of the Rossby 
number. It is therefore of first interest to quantify the 
energy contained in the 2D and 3D components of the 
flow. We estimate the vertically-averaged 2D flow as the 
average of the velocity field over the vertical extent A 2 ; 
of the PIV held, 

1 

U 2 D = ^ / u{x,z)dz, (2) 

Az Jq 

and the remaining z-dependent 3D how as U 3 £i = u — 
U 2 D- We compute the energy per unit mass of these two 
how components as {u^jj)/2 and (u 3 ^)/ 2 , with ~ the 
temporal average and (•) the spatial average over the 
PIV held. They are plotted in Fig. ^a.) as a function 
of the rotation rate fl. Because of the limited height of 
the PIV held, the 2D how estimated from Eq. Q un¬ 
avoidably contains 3D huctuations associated to vertical 
scales larger than Az, so the measured 2D energy may 
overestimate the true one. 

Figurel^a) shows that without rotation the 2D and 3D 
components of the how have comparable energy. With ro¬ 
tation, the 2D energy increases with D, following approx¬ 
imately the power law [4l[, whereas the 3D energy 
remains approximately constant, and represents only 5% 
of the total energy at the largest rotation rate. Although 
most of the energy is contained in the 2D how compo¬ 
nent for D 7 ^ 0 , a signihcant fraction of the dissipation 
still arises from the 3D huctuations. Assuming axisym- 
metry, we compute an estimate of the energy dissipation 
rate e = iy({dui/dxj)^) from the 6 terms of the velocity 
gradient tensor accessible in the 3C2D measurements. 



This dissipation rate, computed both for U 2 d and U 3 D, 
is shown in Fig. [5t^b). Since the derivatives are obtained 
from hnite differences at the smallest resolved scale, the 
computed dissipation underestimates the true one (the 
PIV resolution is 1.75 mm while the Kolmogorov scale is 
of order of 0.6 mm [dlj). However, we expect the mea¬ 
sured evolution of e with D to rehect the true one. 

We hrst compare these bulk energy dissipation rates 
to an estimate of |rictional losses d^e to laminar Ekman 
layers, CEk — v where C/_Lrms is 

°Ek ^ ^ ^ 

the root-mean-square horizontal velocity. This estimate 
ranges from 8 x 10 “® m^ s“^ for D = 2 rpm to 1 x 
10 “^ m^ s“^ for D = 16 rpm: it is smaller than the bulk 
energy dissipation of both the 2D and 3D parts of the 
turbulent flow, by a factor of 10 for slow rotation and 4 for 
rapid rotation. A detailed experimental characterization 
of these Ekman layers would be necessary to validate the 
assumption of laminar layers, but it is beyond the scope 
of the present study. 



FIG. 3. (Color online) Temporal energy spectrum E{a) as 
a function of the angular frequency cr for different rotation 
rates D. (tq indicates the frequency of the opening-and-closing 
cycle of the flaps. The Coriolis frequency 2D is highlighted 
with filled symbols. 

We now compare the bulk energy dissipation rates in 
the 2D and 3D parts of the turbulent flow. Remarkably, 
while the 3D fluctuations represent a small fraction of 
the total energy, they account for a large fraction of the 
dissipation at all rotation rates. It is therefore of first 
interest to investigate these 3D modes, and to determine 
to what extent they can be described in terms of inertial 
waves. 

IV. TEMPORAL ANALYSIS 

We now focus on the temporal dynamics of the velocity 
field, which we characterize through the energy distribu¬ 
tion of turbulent fluctuations as a function of angular 
frequency a. This temporal energy spectrum is defined 
as 

A 'jr 

where 

1 

Ui(x,CT) = — / Mj(x,f)e (5) 

27r Jg 

is the temporal Fourier transform of the velocity held 
Mi(x, t) (with i = x,y,z), T the acquisition duration and 
(•) the spatial average. The normalization is such that 
(rtf) = /p°° E{a) da, with ~ the temporal average. We use 
the standard Welch’s method [i^ to improve the statis¬ 
tical convergence of the power spectrum. 

We plot E{a) for each rotation rate D in Fig. |31 For 
D 7 ^ 0 , we observe a global increase with D of the en¬ 
ergy at all frequencies, consistently with the behavior of 
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FIG. 4. (Color online) Temporal energy spectra of the verti¬ 
cal, i?|| (light gray, blue), and horizoirtal, E±/2 (black), veloc¬ 
ity components for different rotation rates fl. For visibility, 
there is a vertical shift by a factor of 10^ between couples 
of curves at differeirt The dashed lines show power laws 

the overall energy in Fig. [^Ka). These spectra for 17 0 

strongly differ from the non-rotating spectrum, with rela¬ 
tively much more energy at low frequency in the rotating 
case: global rotation induces slow dynamics. 

A first step towards a description of the flow anisotropy 
in the frequency domain can be provided by further de¬ 
composing the power spectrum density (|4]) as 

E{a) = E\\{(t) + EA_{a), ( 6 ) 

with A||((t) = 47 r(|{(z(x, (t)P)/T the spectrum of the 
vertical velocity and E^(a) = 47r(|{ia;(x, cr)p)/T -b 

47 r(|'Uy(x, ct)P)/T the spectrum of the horizontal velocity. 
This decomposition highlights the frequency-dependent 
componentality of turbulence, i.e., the distribution of en¬ 
ergy among the different velocity components, which is 
related to the polarization anisotropy [1^ [TsI . fl^ . This 
is not to be confused with the frequency-dependent di¬ 
mensionality of turbulence, which compares the vertical 



FIG. 5. Componential anisotropy ratio as a function of the 
normalized frequency a* = (t/ 2I7. Isotropy is indicated by 
the horizontal dashed line at 2E\\/E± = 1. The dotted 
line indicates the prediction for a plane inertial wave, i.e., 
2F;||/ifj. = 2(l-a2)/(l + aJ). 


and horizontal characteristic scales at a given frequency 
(described in section IVT). 

The temporal spectra i?||(o') and E±{a)j2 are shown 
in Fig. |4] for all rotation rates. Without rotation, en¬ 
ergy is nearly equally distributed among the 3 velocity 
components (i.e., E± ~ 2ill||). There is actually a slight 
over-representation of horizontal energy, a consequence of 
the forcing device geometry which preferentially injects 
energy in horizontal motions. As the rotation rate in¬ 
creases, the high frequencies remain nearly isotropic (iso¬ 
component), whereas the low frequencies become gradu¬ 
ally anisotropic, with i?||(cr) nearly flat and E±(a) ap¬ 
proaching a power law close to This anisotropy is 

related to the fact that, as 17 increases, the decorrelation 
frequency (i.e., the frequency below which the spectrum 
becomes flat) gets significantly smaller for the horizontal 
velocity (ddecX = 0.04 ±0.01) than for the vertical veloc¬ 
ity (cTdeclb increasing from 0.05 to 0.40 rad s“^ for 17 from 
2 to 16 rpm). For the largest rotation rate (17 = 16 rpm), 
there is a clear range of frequencies over which the hori¬ 
zontal spectrum E±{a) follows a (T“®/^-power-law. This 
range gets narrower for decreasing rotation rate 17. 

The decorrelation frequency of the vertical velocity ap¬ 
pears to scale as CTdecH* = o'dec||/2f7 = 0.25 ± 0.05, which 
becomes evident when plotting the ratio 2E\\/E± as a 
function of the normalized frequency cr* = (7/217 (Fig.[5|). 
In this figure, for cr* > 1, for which no inertial waves can 
exist, energy is nearly equally distributed among the ve¬ 
locity components {2E^^/E± ~ 1). The frequency range 

over which a power law cr* is approached is bounded by 
the two decorrelation frequencies: on the left by crdecj. 
and on the right by crdec|| ■ Interestingly, we also ob¬ 
serve a small frequency domain crdec||* < cr* < 0.6 over 
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which energy in the vertical component is slightly larger 
than in each horizontal component. Such a slight over¬ 
representation of the vertical velocity is compatible with 
the “componential” anisotropy of an assembly of linear 
inertial waves: a single plane inertial wave has a compo¬ 
nential anisotropy 2E\\/Ea_ = 2(1 — crj)/(l -I- cr*) (shown 
as a dotted line in Fig. ED, with a larger rms velocity 
along the vertical than along any horizontal direction for 
a* < 1/V3~ 0.6. 

As can be seen in Fig. [21 for the experiment under 
rapid rotation, nearly 90% of the energy is contained in 
the 2D vortex mode. As discussed quantitatively in Cam- 
pagne et al. [2^ , this strong 2D nature of the flow drives 
an inverse cascade of energy for the horizontal velocity 
and a direct cascade of energy for the vertical velocity. 
The horizontal velocity consequently exhibits slow dy¬ 
namics, while the vertical velocity fluctuations are found 
at higher frequencies (Fig. |4D. This behavior is consis¬ 
tent with the usual phenomenology of rapidly rotating 
turbulence: the flow becomes approximately 2D at low- 
frequency, and the vertical velocity behaves as a passive 
scalar, which is stretched and folded by the horizontal 
velocity. This produces thin vertical layers swept by the 
horizontal flow, yielding rapidly changing time series of 
the vertical velocity. In a similar fashion, the 
power-law of E±{a) could originate from the stochastic 
sweeping by the large scale horizontal flow of a spa¬ 

tial spectrum, reminiscent of the inverse energy cascade 
of 2D turbulence. 


V. SPATIO-TEMPORAL ANALYSIS 
A. Spatio-temporal correlations 

We now turn to a combined spatio-temporal analysis of 
the PIV time series, focusing on the signature of inertial 
waves in terms of dimensional anisotropy. This signature 
is sought here in terms of characteristic horizontal and 
vertical scales of the turbulent structures as a function 
of their frequency. For this, we define the frequency- 
dependent two-point spatial correlation of the temporal 
Fourier transform of the velocity field 

OiTT 

^) = Y + r, ct) + c.c.), (7) 

with * the complex conjugate (here again Welch’s method 
is used to improve convergence). Instead of the spectra 
considered in Refs. [H, we compute spatial correla¬ 
tions, because the former are sensitive to finite size effects 
arising from the PIV field being of limited extent com¬ 
pared to the largest flow structures. The correlation © 
probes the energy distribution among vector separations 
r for each frequency a. Its angular average provides an 
estimate for the cumulative energy from scale r = |r| 
to r = oo for turbulent motions of frequency cr. The 
single-point limit of this correlation is the temporal en¬ 
ergy spectrum Q, i.e. E{a) = R{r = 0,cr). 



0 30 60 90 120 0 30 60 90 120 
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FIG. 6. (Color online) Maps of the normalized two-point 
correlation R(r,a)/E{a) in the vertical plane (rx,r||) for 
D = 0 (left) and D = 16 rpm (right), at four frequen¬ 
cies a = 0.1,0.7, 2.4, 3.3 rad s“^. In the rotating case, 
the corresponding normalized frequencies are ct* = a/2Q. = 
0.03,0.20,0.71,0.98. Isocontour lines for R{r,a)/E{a) < 0.2 
are noisy and are not shown. 
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Maps of the normalized correlation R{r,a)/E{a) are 
plotted in Fig. [5] for = 0 and 16 rpm at four selected 
angular frequencies a. In the non-rotating case, the iso-i? 
lines are approximately circular at all scales and frequen¬ 
cies, indicating the overall isotropy of turbulence. The 
strongly peaked correlation that develops around r = 0 as 
a increases indicates that rapid turbulent fluctuations are 
found at small scales only. In the rotating case, the iso- 
R lines evolve gradually from quasi-vertical at small fre¬ 
quency (“cigar” anisotropy) to more horizontal for cr* ^ 1 
(“pancake” anisotropy), with cr* = aj2n the normalized 
frequency. The “cigar” anisotropy observed at cr* <C 1 is 
consistent with the 2D3C vertical invariance predicted by 
the Taylor-Proudman theorem for vanishing frequency: 
it corresponds to the zero-frequency limit of Eq. ([T]) for 
nearly horizontal wave vector. Similarly, the tendency 
towards “pancake” anisotropy, observed for cr* —>■ 1 and 
sufficiently large r ±, is consistent with the nearly vertical 
wave vector limit of Eq. O- 

A natural way to characterize the frequency-dependent 
anisotropy would be to compute integral scales along and 
normal to the rotation axis at each frequency. Here we 
consider a finer approach, which also takes into account 
the scale-dependence of this anisotropy: for each fre¬ 
quency (T and horizontal scale r±, we identify the ver¬ 
tical scale •^||(rx,o') at which the correlation along the 
vertical axis is equal to the one at r±e±, i.e. such that 
i?(r = ^||e||,tT) = R{r = r±e±,a). In practice, we com¬ 
pute as the vertical semi-axis obtained from the fit of 
the iso-i? line defined by i?(r, a) = i?(rxex, cr) with an el¬ 
lipse of prescribed horizontal semi-axis r±_. This method 
allows us to filter out the noise in the iso-i? lines at small 
R. It also extends the analysis to values of larger than 
the PIV field height (Az = 140 mm), which is useful at 
small (T* for nearly vertically invariant i?. Poor fits de¬ 
fined by a correlation coefficient less than 0.9 or such that 
i|| is larger than 2Az are discarded. We finally define the 
scale and frequency dependent anisotropy factor as 


A(rx,cr) 


^||(»'-L,cr)’ 


( 8 ) 


It is equal to 1 for isotropic turbulence, to 0 for vertically 
invariant (2D3C) turbulence, and to oo for horizontally 
invariant (1D2C) turbulence. 

If the anisotropy of the two-point correlation i? at fre¬ 
quency cr* ^ 1 is governed by linear inertial waves, we 
expect A to be independent of the scale, and to be set by 
the dispersion relation O. A simple estimate, assnming 
that the wave vector k is in the vertical measurement 
plane and identifying on dimensional grounds ^ 

and r± ^ kj^, yields a frequency-dependent anisotropy 
factor — 1)“^^^. Considering now an as¬ 

sembly of inertial waves with an axisymmetric wavevec- 
tor distribution, an analytic computation (given in the 
Appendix) leads to a very similar result. 


Aiw{<^*) 



(9) 
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FIG. 7. Anisotropy factor A ([8]) as a function of the normal¬ 
ized frequency at — (t/ 2H at different rotation rates (same 
symbols as in Fig. El), for three horizontal scales r± = 9, 18 
and 50 mm. The continuous line represents the inviscid 
inertial-wave prediction Aiw 


In Figs. [7](a-c), we compare the anisotropy factor A 
measured in the four rotating experiments to the inertial- 
wave prediction m at three horizontal scales r± = 9, 18 
and 50 mm. For all scales and cr* < 1, A is an increas¬ 
ing function of cr*, confirming that slow fluctuations are 
more vertically elongated than fast fluctuations. We find 
that the inertial-wave prediction (O provides a good de- 
















scription of the data at large horizontal scales and large 
rotation rate. For such large scales {r± ~ 50 mm), the 
anisotropy factor is no longer accessible for cr* <0.1 be¬ 
cause it corresponds to much larger than the height of 
the PIV field. On the other hand, at smaller horizontal 
scale the prediction ([9|) fails, with small frequencies more 
isotropic than predicted by the inertial-wave argument. 

Because of the moderate Reynolds and Rossby num¬ 
bers of the present experiment, two effects may be consid¬ 
ered to explain why large scales follow the inertial-wave 
prediction whereas small scales do not: viscous damp¬ 
ing and sweeping of small scales by the velocity at larger 
scales. Viscosity introduces an imaginary term 
in the dispersion relation o without modifying its real 
part. Waves (cr, k) such that |k| ^ r~^ are therefore 
damped, with r^, = \fvja a viscous cutoff. This vis¬ 
cous cutoff is of order of 1 to 10 mm for the normalized 
frequencies cr* = CT/2f2 in the range [10“^,1] considered 
in Fig. [5l However, since viscous damping affects the 
wave amplitude without modifying the wave vector com¬ 
ponents, its should not affect the anisotropy. We there¬ 
fore focus in the following on the sweeping effect. 

B. Sweeping effect 



FIG. 8. Rescaled intrinsic frequency V m as a function of 
the sweeping parameter S. The symbols indicate the different 
rotation rates, and are the same as in Fig. [S] The dashed 
line N 1/S shows the low-S prediction for (non-swept) 
ensembles of inertial waves [Eq. (I12II ]. For S <C 1 one sees the 
signature of the dispersion relation o, while for S' S> 1 one 
detects swept inertial waves. 


Sweeping corresponds to the advection of the waves by 
the large-scale flow, which leads to a modification of their 
apparent frequency. An inertial wave propagating in a 
time-independent uniform flow U has a Doppler-shifted 
frequency. 


cr = (Ti + k • U , 


( 10 ) 


where ai is the intrinsic frequency given by o, and a 
is the frequency at which the wave is detected in the 
frame of the rotating tank. In our experiment, the en¬ 
ergetic large-scale 2D flow may be thought of locally as 
a uniform sweeping flow U that evolves slowly in time, 
inducing a scrambling of the waves’ spatio-temporal sig¬ 
nature. An order of magnitude of the typical Doppler- 
shift can be estimated by k±U±rms, where U±rms is the 
root-mean-square horizontal velocity. A key difference 
between equations m and © is that the frequency a 
now depends on the magnitude of k, with small-scale 
waves more affected by sweeping. 

For an ensemble of inertial waves with axisymmetric 
wavevector statistics, the intrinsic frequency ai can be 
related to the anisotropy through Eq. ®. Subsituting 
the corresponding expression into (US and estimating the 
Doppler-shift term on dimensional grounds, we obtain 


2D 


Vl + 2A-2 


+ C 


U 


Lrms 

r± 


( 11 ) 


where C is a constant of order unity. This indicates that 
the parameter 


N = 


2nr± 



( 12 ) 


should be a unique function of the sweeping parameter 
S = U±rms/o'r±. N corresponds approximately to the 
intrinsic frequency of inertial waves rescaled by the ad- 
vective time r±/U±rms, whereas S is the observed period 
of the waves, rescaled by the same advective time. 

Figure [5] confirms this picture: the data for differ¬ 
ent values of D, r± and a collapse onto a master curve 
N = f{S). This collapse indicates that sweeping is in¬ 
deed responsible for the departure from the inertial wave 
prediction at small frequencies and/or small scales. The 
expected asymptotic behavior for small sweeping param¬ 
eter is ~ 1/S, which corresponds to the prediction 
(P|) for an axisymmetric ensemble of non-swept inertial 
waves. The data is in quantitative agreement with this 
small-5' prediction, shown as a dashed line in figure[8l For 
large 5, equation (HU indicates that N should asymptote 
to a constant value N ~ C, which again is compatible 
with the data. 

The master curve in figure [5] has the following simple 
interpretation: high-frequency or large-scale waves are 
hardly affected by sweeping. The Doppler-shift term is 
negligible compared to their intrinsic frequency, and their 
location in a space-time energy distribution is given by 
the dispersion relation ©■ This is the low-5 behavior in 
figure [5] 

By contrast, when focusing on low frequencies a or 
small-scales in the frame of the tank, one measures the in¬ 
ertial waves with intrinsic frequency at = a, but one also 
detects many waves with at ^ a that are Doppler-shifted 
back to frequency a by the advective term in (ITO)) . The 
anisotropy measured at low-frequency a therefore results 
from strongly swept inertial waves with various intrin¬ 
sic frequencies, and the information from the dispersion 
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relation o is lost in the space-time correlation. The 
limit S' ^ 1 corresponds to frequencies a that are much 
lower than the inverse advective time. In this cr —> 0 
limit, one detects mostly waves with Ui ^ a that are 
Doppler-shifted by the horizontal flow in such a way that 
they are almost steady in the frame of the tank: this is 
a cr-independent regime that corresponds to the large-S 
plateau in figure [51 


Aubertin, L. Auffray, C. Borget and R. Pidoux for their 
experimental help. This work is supported by the ANR 
grant no. 2011-BS04-006-01 “ONLITUR”. B.G. acknowl¬ 
edges support from Labex PALM. F.M. acknowledges the 
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Appendix A: Anisotropy factor for a statistically 
axisymmetric distribution of inertial waves 


VI. CONCLUSION 

In the present experiment, the anisotropy of the turbu¬ 
lent energy distribution at a given spatio-temporal scale 
(rx, cr) is well-described by the inertial-wave dispersion 
relation at high-frequency and/or large-scale only. The 
smaller-scale waves are subject to intense sweeping by 
large-scale turbulent motions contained predominantly 
in the 2D “vortex” mode. This conclusion is compatible 
with the numerical findings of Clark di Leoni et al. [s^, 
who also identify the sweeping time scale as the relevant 
decorrelation time at small scale. 

Such sweeping by the 2D mode has strong implica¬ 
tions for wave-turbulence theories. Indeed, most waves 
do not follow the inertial-wave dispersion relation, and 
the assumptions of weak turbulence theory break down 
even at the linear stage in wave amplitude: instead of 
the dispersion relation the linear problem consists in 
determining the evolution of waves embedded in a tur¬ 
bulent 2D flow. This is a formidable task in general, 
because the 2D flow is space- and time-dependent: in 
the discussion of our data, we simplified the problem by 
assuming that the 2D flow is at much larger scales and 
slower frequencies than the waves, therefore including it 
as a simple Doppler-shift term in the dispersion relation. 

We conclude with a discussion on the dimensionality of 
the forcing. In the present experiment, the flow is driven 
by vertically invariant flaps: such a quasi-2D forcing de¬ 
vice enhances two-dimensionalization and the resulting 
sweeping of the 3D flow structures. Nevertheless, accu¬ 
mulation of energy in the 2D mode is a robust feature of 
rotating turbulence, that takes place for arbitrary forc¬ 
ing geometry, even if the forcing does not input energy 
directly into the 2D mode. A careful and extensive nu¬ 
merical study of this issue has been recently reported for 
the fully-3D Taylor-Green forcing [d^: for rapid global 
rotation and low viscosity, energy accumulates in the 2D 
mode until the Rossby number based on the velocity of 
this 2D flow is of order unity. If these findings are con¬ 
firmed, the sweeping of the most energetic 3D structures 
would be an inevitable outcome of this accumulation of 
2D energy. 
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We compute the anisotropy factor A for an ensemble 
of independent plane inertial waves, with axisymmetric 
wave vector statistics. The temporal Fourier transform 
of the velocity field reads 

u(x, 0 -)= f a(k, cr)e**''’^dk, (Al) 


where a(k, cr) is the space-time Fourier amplitude of the 
velocity field at wave number k and frequency a. The 
two-point velocity correlation at frequency cr © can be 
written 

^(r^cr) = i ^JJa(ki,cr) • a*(k 2 ,cr) 

X e*(‘^i'^”‘^^'(^+'’»dkidk 2 -kc.c.^ , 

= J |a(k, ct)P cos(k • r)dk, (A2) 

where (•) is the space average, and r is a separation vec¬ 
tor inside the PIV plane. Introducing spherical coordi¬ 
nates with vertical polar axis, we denote ip the azimuthal 
angle between k and the vertical PIV plane. The argu¬ 
ment of the cosine becomes 


k • r = fc|| r|| -I- k±r± cos (p . (A3) 


Let us first consider an ensemble of inertial waves hav¬ 
ing the same wavenumber, |k| = k. For a given reduced 
frequency cr* = aj2n, the dispersion relation determines 
the ratio |A:|| |/fcx, and because of statistical axisyinmetry 
|a(k, cr)p is independent of ip. The spatial correlation at 
frequency a becomes 


R{v,a)= 


G{k,c 

27r 


r2Tr 


cos(A;||r|| + k±r± cos ip)dip, 


G{k, cr) cos (cr*fcr||) Jq ^\/1 — crj kr±^ , (A4) 


where Jq is the Bessel function of the first kind and 
G{k, cr) is a prefactor proportional to the squared am¬ 
plitude of the waves at wavenumber k and frequency cr. 

For a given frequency cr and horizontal scale r±, the 
vertical scale ^\\(r^,a) is determined from the isolines 
R =constant in the (rxiTu) plane. An isoline of R start¬ 
ing on the horizontal axis at r± intersects the vertical 
axis at r|| = £||(rx,cr). According to expression (IA4I) . 
such isolines connecting the two axes exist provided the 
argument of the Bessel function is smaller than the first 
zero of this function, i.e.. 


\/l - crj kr± < Go , 


(A5) 
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where Jo(C'o) = 0 (Cq — 2.40). Equating expression (jA4l) 
computed for (r_L,r|| = 0) and for (rj_ = Ojrn = ^y) leads 
to 


arccos 


and an anisotropy factor 
A(r_L,cr*) = - 


Jo (\/ 


1 - 0-2 


a^k 

kr±a:^ 


arccos 


Jo (y 


1 — 0-2 


(A6) 


(A7) 


This anisotropy factor depends very weakly on kr±: it is 
minimum for low kr±, where Taylor expansion for kr± <C 
1 leads to 


A(r_L, o-*) ~ -\/2 



(A8) 


whereas it is maximum for the maximum value of r± 
allowed by (ESI), which gives 




2Cq (7^ 

^ \/r^ 


1.53 


CT ^ 

\/l - 


(A9) 


Because the numerical prefactors in (IA8I1 and (IA9I1 differ 
by less than 10%, we can say that the anisotropy factor 
of this statistically axisymmetric distribution of inertial 
waves is given by expression (IA8I) within 10% accuracy. 
Because this anisotropy factor is almost independent of fc, 
we finally expect it to be approximately given by expres¬ 
sion (IA8I) for a realistic superposition of inertial waves 
with different wave numbers k. 
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